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ABSTRACT 

Recently Sahu et al., using the Hubble Space Telescope to monitor stars in the direction of the old 
(~ 12 Gyr) globular cluster M22, detected six events in which otherwise constant stars brightened by 
~ 50% during a time of ^ 1 day. They tentatively interpret these unresolved events as due to microlensing 
of background bulge stars by free-floating planets in M22. Using simple analytic arguments, I show that 
if these spike events are due to microlensing, the lensing objects are unlikely to be associated with M22, 
and unlikely to be part of a smoothly distributed Galactic population. Thus either there happens to 
be a massive, dark cluster of planets along our line-of-sight to M22, or the spike events are not due 
to microlensing. The lensing planets cannot be bound to stars in the core of M22: if they were closer 
than ~ 8 AU, the lensing influence of the parent star would have been detectable. Moreover, in the 
core of M22, all planets with separations > 1 AU would have been ionized by random stellar encounters. 
Most unbound planets would have escaped the core via evaporation which preferentially affects such 
low-mass objects. Bound or free-floating planets can exist in the outer halo of M22; however, for 
reasonable assumptions, the maximum optical depth to such a population falls short of the observed 
optical depth, t~3x 10~ 6 , by a factor of 5-10. Therefore, if real, these events represent the detection of 
a significant free-floating Galactic planet population. The optical depth to these planets is comparable 
to and mutually exclusive from the optical depth to resolved events measured by microlensing survey 
collaborations toward the bulge, and thus implies a similar additional mass of lensing objects. If the 
spatial and kinematic distributions of the two populations are the same, there are > 10 3 planets per 
bulge microlens. Such a population is difficult to reconcile with both theory and observations. 

Subject headings: planetary systems — globular clusters: individual (M22) — gravitational lensing 



1. INTRODUCTION 

Our understanding of the low-mass end of the compact 
object mass function has seen considerable progress on 
many fronts in the last decade. The enormous success 
of radial velocity searches for extrasolar planets has im- 
proved considerably our knowledge of the statistics, dis- 
tribution, and mass function of close (a <^ 3 AU) and 
relatively massive (M ;> 0.2 Mj) companions to nearby 
stars (Marcy, Cochran, & Mayor 2000; Jorissen, Mayor, & 
Udry 2001; Tabachnik & Tremaine 2001). Direct searches 
for brown dwarf companions to local solar-type stars have 
led to only one detection, GL 229B (Nakajima et al. 1995, 
Oppenheimer et al. 2001). However, serendipitous discov- 
eries in wide-field surveys of very low mass companions to 
normal stars (Kirkpatrick et al. 2001) has led to the con- 
clusion that the 'brown dwarf desert,' the paucity of close 
brown dwarf companions to stars monitored in radial ve- 
locity surveys (Marcy & Butler 2000), does not exist for 
wide separations (Gizis et al. 2001). The statistics of very 
low mass companions M <; 10 M®, will likely have to wait 
for future astrometric (Lattanzi et al. 2000) or transit- 
detecting (Borucki et al. 1997) satellites. 

These studies have focused on GKM dwarfs in the im- 
mediate local neighborhood. Detailed studies of the parent 
stars of extrasolar planets have revealed that these hosts 
have higher metallicity in comparison to an unbiased field 

bubble fellow 



sample (Gonzalez 1997; Laughlin 2000). However, it is 
difficult to interpret this observation, as it is not clear if 
the cause of this enhanced metallicity is stellar pollution 
from cannibalized planets or rather that low metallicity 
tends to prohibit formation (Murray et al. 2001; Santos, 
Israelian, & Mayor 2001; Pinsonneault, DePoy, & Coffee, 
2001). In order to resolve this issue, it may be necessary 
to look toward other systems, such as globular clusters, 
open clusters, and the Galactic bulge. Such systems are 
useful in that they act as a control samples, where due 
to the homogeneous nature of the systems, one particular 
parameter that may affect planet formation and/or evo- 
lution (i.e. metallicity, age, local density) can be better 
isolated. One of the most promising methods of studying 
companions to stars in such systems is transits, which is 
in a sense ideally suited to this application since it works 
best when a large number of stars can be monitored simul- 
taneously. An important null result was found when this 
tactic was combined with the high resolution of the Hubble 
Space Telescope (HST) to search for planetary compan- 
ions of stars in the (relatively) metal-rich globular cluster 
47 Tuc. Gilliland et al. (2001) found that the frequency 
of close planetary companions to stars in 47 Tuc is more 
than an order of magnitude smaller than in the local so- 
lar neighborhood. This lack of planets is difficult to ex- 
plain in terms of dynamical stripping (Davies & Sigurdsson 
2001) but may be explicable by invoking disk photoevap- 
oration (Armitage 2000). Observations in other environ- 
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mcnts, such as the Galactic bulge (Gaudi 2000) or open 
clusters, have been considered, are currently being under- 
taken, or are being planned. 

Large scale optical and infrared surveys such as the 
Sloan Digital Sky Survey (SDSS) and the Two Micron All 
Sky Survey (2MASS) have discovered a significant popu- 
lation of low mass objects. By now nearly 100 L-dwarfs 
(Fan ct al. 2000, Kirkpatrick et al. 2000), and a dozen 
T-dwarfs (Strauss et al. 1999, Burgasser et al. 1999) are 
known. The difficulty with using these observations to 
constrain the low mass field object mass function is that 
for a given effective temperature, there exists a degeneracy 
between mass and age, making the determination of the 
mass function fundamentally uncertain (Reid et al. 1999). 

Observations of low mass field objects confined to stellar 
aggregations alleviate some of these difficulties, because 
the age of the system is known or can be estimated in 
many cases. Such studies have led to the detection of 
a population of low-mass, free-floating objects in several 
open clusters (Bouvier et al. 1998; Hambly et al. 1999; Za- 
patero Osorio et al. 2000; Martin et al. 2000). Many of 
these are clearly brown dwarfs, but some are suggested to 
have masses that place them in the canonical planet regime 
(M < 13Mj). The relation of these objects to those found 
orbiting local stars is unknown, and even their designation 
as planets is controversial (McCaughrean et al. 2001). 

Microlensing is a unique method of detecting planets 
that offers the significant advantage that its effect is sen- 
sitive only to the mass of the lensing object, and does 
not rely on either the flux of the planet or its parent star. 
Therefore distant, extremely faint, or even completely non- 
luminous compact objects, either free-floating or bound 
to other objects, can be detected using microlensing. In- 
deed microlensing was originally suggested by Pacyhski 
(1986) as a method to look for dark matter in the halo of 
our Galaxy. In principle, extremely low-mass objects can 
be detected, with the ultimate lower limit set by the size 
of the source stars. This fact was exploited by the MA- 
CHO and EROS collaborations to rule out objects having 
masses of 10~ 7 — 10 _3 M Q as the primary constituents of 
halo dark matter (Alcock et al. 1998). The PLANET col- 
laboration acquired and analyzed five years of photometric 
data searching for planetary companions to Galactic bulge 
stars, using a method first suggested by Mao & Paczyhski 
(1991), but found no candidates, implying that ~ Mj mass 
planets with separations of 1.5 AU < a < 4 AU occur 
in less than 1/3 of systems (Albrow et al. 2001, Gaudi 
et al. 2002). Free-floating planets are also detectable in 
microlensing searches toward the Galactic bulge (di Ste- 
fano & Scalzo 1999a). The primary drawback to using 
microlensing to detect low-mass objects is that the mass 
of the lensing object typically cannot be directly measured. 
Instead, one measures the timescale £e, which is a degen- 
erate combination of the mass of the lens, and the lens- 
source relative parallax and proper motion. This makes 
unambiguous detection of low-mass objects difficult. 

Recently, Sahu et al. (2001, hereafter SM22) presented 
an ingenious method of overcoming this difficulty, based on 
a suggestion originally made by Paczyhski (1994). They 
used HST to monitor stars toward the old, metal-poor 
globular cluster M22 (I = 9°.9, b = -7°.6). Because M22 
is projected in front of the bulge, there exists a significant 



probability (optical depth) that background bulge stars 
will be microlensed by foreground M22 objects. Since the 
distance, velocity dispersion, and proper motion of M22 
are known, the detection of a microlensing event, and a 
measurement of £e yields a measurement of the mass of 
the lens, with only the unknown distance and velocity of 
the source contributing significantly to the error. In fact, 
SM22 did detect one event with £e = 17.6 days, yielding 
a lens mass of M — O.IS^q'^Mq - an impressive mass 
measurement of an M-dwarf 2.6 kpc away. 

SM22 also detected six events in which a constant light 
curve brightened by ~ 50% during one set of two mea- 
surements, separated by six minutes, and then returned 
to baseline for the remaining measurements. By direct 
inspection of the images, SM22 rule out detector arti- 
facts or cosmic rays as the source of these brightenings. 
Furthermore, the brightening was consistent in both im- 
ages in each of the six cases. SM22 rule out several other 
astrophysical sources of variability, and tentatively con- 
clude that these brightenings are due to unresolved mi- 
crolensing events. With the sampling rate of ~ 1 per 
4 days, this implies an upper limit to the timescale of 
t~E 1 day. If the lenses are located at the distance of M22 
[Di = 2.6±0.3 kpc (Peterson & Cudworth 1994)] and have 
kinematics similar to M22 stars, then this translates to an 
upper limit on the mass of the lenses of M p <; 0.25Mj. 
SM22 tentatively conclude that they have found evidence 
for a considerable population of free-floating planets in 
M22. This is interesting: if every star in M22 had n p plan- 
ets of mass M p associated with it, than the optical depth 
t p to lensing by these planets would be roughly given by 

where M* is the mass of a typical star in M22, and t* is 
the lensing optical depth to normal stars. Therefore for 
M p < 0.25Mj and t* ~ 10~ 5 for the core of M22 (see 
§5), T p <; 10~ 8 n p . Since SM22 monitored ~ 8 x 10 5 stars, 
detecting any events would be very unlikely unless n p is 
very large. Indeed, SM22 estimate that 10% of the cluster 
mass may be in these planetary objects. 

de la Fuente Marcos & de la Fuente Marcos (2001) dis- 
cuss some of the implications of this inference. By ana- 
lyzing a suite of N-body simulations, they argue that the 
majority of planets in the core of M22 will be ionized from 
their parent stars, and those ionized planets will quickly 
evaporate from the core. They argue that bound planets 
may be able to survive in the outskirts of M22, thus repro- 
ducing the observed optical depth, but only if multiplanet 
systems are common. However, they favor the possibility 
that these events are due to a dark cluster of planets not 
associated with the globular cluster. 

Because microlensing events do not repeat, the only way 
to directly verify these observations is to perform the same 
experiment at even higher time resolution. Since this ex- 
periment essentially requires the resolution of HST, this 
would be a major investment of limited resources. Given 
the importance of this result on the one hand, and the 
considerable expenditure of resources required for direct 
verification on the other, it seems prudent to consider the 
SM22 interpretation of these spike events in more detail. 
Here I present such a study, which is complementary to 
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sources in the literature that will be relevant to the dis- 
cussion. In some cases, there are discrepancies in the pub- 
lished literature. Where possible, I have adopted the most 
modern determinations. Adopting other values for these 
parameters will change some of the resulting computations 
in detail, but will not affect the general arguments and 
conclusions. 

3. MICROLENSING BASICS 

In this section I review the fundamentals of microlens- 
ing, concentrating on those concepts that are important 
for the present study. A lens of mass M at a distance Di 
has an angular Einstein ring radius of 
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and D s is the distance to the source. The scaling relation 
at the extreme right hand side of equation (^), and all scal- 
ing relations presented in this section, are appropriate for 
M22 parameters, i.e. Dj = 2.6 kpc, D t = 8.0 kpc. At the 
distance of the lens, #e corresponds to a physical distance 
of 
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that presented by de la Fuente Marcos & de la Fuente 
Marcos (2001) in that I rely primarily on analytic argu- 
ments (rather than on the results of N-body simulations). 
In §2, I briefly review the SM22 observations, and collect 
some relevant parameters for M22. In §3, I review the 
basics of microlensing, and apply these to the M22 ob- 
servations in §3.1. In §4, I show that these objects are 
unlikely to be associated with M22, by first demonstrat- 
ing that these planets must be separated by a ;> 8AU from 
their parent stars (§4.1), by second showing that planets 
in the core with separations a > 0.3AU would be ion- 
ized by other stars (§4.2), by third showing that most low 
mass, free-floating planets would have evaporated from the 
cluster core over the cluster lifetime (§4.3), and finally by 
showing that the maximum optical depth to planets in the 
halo of M22 falls considerably short of the observed opti- 
cal depth (§4.4). Thus if the spike events are indeed due 
to microlensing, the lensing objects cannot be associated 
with M22. In §5, I explore the implications of this conclu- 
sion, and show that, if smoothly distributed, the implied 
free-floating Galactic planet population is hard to recon- 
cile with current theoretical and observational constraints. 
Thus either the explanation of de la Fuente Marcos & de 
la Fuente Marcos (2001) is correct - there is a dark cluster 
of planets along our line-of-sight to M22 - or the major- 
ity of the spike events cannot be due to microlensing. In 
light of the circumstantial evidence against the microlens- 
ing interpretation of these events, in §6 I present a critical 
re-evaluation of the direct evidence in favor of microlens- 
ing. I summarize and conclude in §7. 

2. OBSERVATIONS AND M22 PARAMETERS 

SM22 monitored the central region of M22 with HST 
over the course of 114 days, with a total of 43 epochs with 
three fields monitored at each epoch. The majority of the 
measurements was taken in the F814W (/) filter with every 
fourth measurement in the F606W (wide V) filter. I will 
be primarily considering the light curves presented in Fig- 
ure 2 of SM22, which show the F814W band light curves 
of the six spike events over a period of 105 days, with ~ 25 
measurements in each light curve, for a sampling interval 
of ~ 4 days. The typical photometric errors for these light 
curves, as judged by the scatter, is 1 — 5%. Table 1 gives 
the the baseline F814W and F606W magnitudes of these 
six spike events, reproduced from Table 1 of SM22, and 
retaining their lettering scheme. Also presented in Ta- 
ble 1 are estimates of the amount by which these events 
brightened, both in magnitudes, AmpsiAw and the corre- 
sponding magnification A assuming no blending 2 . Values 
for AmF8i4 have errors of 2 — 7%, as estimated by eye. 
Finally, Table 1 shows the dereddened color and absolute 
magnitude of the SM22 sources, assuming they are located 
at the Galactic center, i.e. a distance modulus of 14.52 
(8 kpc), and E(B — V) = 0.326. The latter was found us- 
ing the reddening maps of Schlegel, Finkbeiner, & Davis 
(1998). This value agrees well with that reported by Pi- 
otto & Zoccali (1999) from main sequence fitting to M22, 
and translates to Afqq6W — 0.94 and Apsuw = 0.63. 

Table 2 lists parameters of M22 collected from various 

2 For typical ground-based resolutions (~ 1"), a large fraction of sources are blended, i.e., there exist multiple stars in each resolution element 
(Han 1997a). However, with the resolution of HST, essentially all stars are resolved (Han 1997b). Thus the observed flux is most likely directly 
related to the magnification, unless the lens itself, or an unlenscd companion to the source or lens is of comparable brightness to the lensed 
source. 
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The timescale of a microlensing event is given by, 
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where ^t ro i is the lens-source relative proper motion. I 
will assume that the mean source proper motion is zero 
(which is appropriate for sources in the bulge), and thus 

Mrei = Mi = 10-9 mas Y r_1 - 

The foreground lens magnifies the background source by 
an amount that depends only on the angular separation 9 
between the lens and source in units of 9e, u = 9/9e- For 
a single lens the magnification is, 
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where to is the time of maximum magnification and uq is 
the impact parameter. Note that uq is distributed uni- 
formly. 

The optical depth to microlensing, r, defined as the 
probability that any star is magnified by > 1.34, is given 

by 
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Source 








A 

Si 


(m,F606W — mF8i4w)o l 


MF8UW 


Fie- 1 


18.86 


17.85 






0.70 


2.70 


A 


23.32 


21.64 


-0.46 


1.53 


1.37 


6.49 


B 


23.81 


22.25 


-0.46 


1.53 


1.25 


7.10 


C 


23.51 


22.12 


-0.72 


1.94 


1.08 


6.97 


D 


22.70 


21.33 


-0.31 


1.33 


1.06 


6.18 


E 


23.03 


21.49 


-0.49 


1.57 


1.23 


6.34 


F 


22.14 


20.81 


-0.31 


1.33 


1.02 


5.66 



Table 1 Measured and estimated parameters for the M22 events. 

1 Dereddened color assuming E[B - V) = 0.634, or E(F&0&W - F&1AW) = 0.31. 



2 Absolute magnitude assuming a distance modulus of 14.52 and E(B — V) = 0.634. 



Parameter 


Symbol 


Value 


Reference 


Distance 


Di 


2.6 kpc 


Peterson & Cudworth (1994) 


Core Radius 


Tc 


l'A 


Trager, King, & Djorgovski (1995) 


Tidal Radius 


n 


28'. 9 


Trager, King, & Djorgovski (1995) 


Central Density 


Po 


4.65 x 10 3 M Q pc" 3 


Peterson & King (1975) 


Velocity Dispersion 


a 


11.4 km s _1 


Peterson & Cudworth (1994) 


Proper Motion 




10.9 mas yr _1 


Peterson & Cudworth (1994) 


Age 


T 


12 Gyr 


Davidge & Harris (1996) 


Typical Mass 


M* 


0.33 M e 


Paresce & De Marchi (2000) 


Reddening 


E{B - V) 


0.326 


Schlegel, Finkbeiner, & Davis (1998) 


Galactic Coordinates 


(l,b) 


(9°9, -7°6) 




Metallicity 


[Fe/H] 


-1.54 ± 0.11 


Lehnert, Bell, & Cohen (1991) 



Table 2 Collected relevant parameters for M22. 



where p is the mass density along the line of sight to the 
source. 

Observationally, the optical depth can be estimated 
from a sample of N e events with measured timescales t^ % 

by 

7T N ° t 

Tabs = ■ ^2 ~> (Resolved Events) (9) 

2 — 1 

where iV* is the number of stars monitored for a duration 
At, and Ci is the detection efficiency of event i. As I is 
show in §3.1, the expression for unresolved events (where 
£b is not known) takes a slightly simpler form. 

3.1. Application to the M22 Events 

Before discussing constraints on the location of the 
lenses, I first apply the microlcnsing formalism just pre- 
sented to the M22 events. SM22 reports that the single 
resolved event has a timescale of t^ = 17.6 days. As- 
suming that the lens is associated with M22, and that 
the source is in the bulge, this translates into a mass of 
M = O.13M [Eq. p|. Careful accounting of the dispersion 
in bulge source distances and kinematics establishes an er- 
ror of - 15% (SM22). Inspecting the mass function (MF) 
of M22 as determined by Paresce & De Marchi (2000), 
the number of stars with this mass is approximately 20% 
smaller than at the peak of the MF (M = O.33M ). One 
would expect - 0.8(0. 13/0.33) 1 / 2 ~ 50% fewer events from 
objects with M = 0.13M Q than the peak of the MF, and 



thus the derived mass is quite plausible. The optical depth 
implied by this one resolved event ultimately depends on 
the detection efficiency, but a lower limit can be found by 
assuming e = 1. Due to the high quality of the HST pho- 
tometry, the fact that the sampling interval is considerably 
less than £e, and that blending is likely not an issue, the 
true optical depth is probably not much larger than this 
lower limit. I find, 

r res £ 3 x 10~ 6 . (10) 

Where I have adopted V* = 83, 000, as indicated by SM22, 
and At = 105 days, which is derived from their Figures 1 
and 2. As I show in §4.4, this lower limit is a factor of 
~ 3 smaller than would be expected based on what we 
think we know about the central density and structure 
of M22. This could be caused by inaccurate parameters 
for M22, an efficiency considerably smaller than unity, or 
Poisson fluctuations. The one resolved event is somewhat 
unusual because the source is reported to be a variable. In 
fact, in order to derive £e from the light curve, SM22 fit 
the event to a binary-source model. This provides a con- 
sistency check on the microlensing interpretation, as the 
binary-source model provides a limit on #e which can be 
compared with the value derived from the timescale (Han 
& Gould 1997, Alcock et al 2001). 

Obviously the timescales of the six unresolved events are 
unknown. SM22 assert that the statistical upper limit to 
tE is 0.8 days (at a 95% confidence level), which translates 
to M < 0.27 x 10~ 4 M Q or M <; 0.25Mj where Mj is the 
mass of Jupiter. A lower limit to the mass can be derived 
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by considering finite source effects, and is M ^> 3 x 10~ 8 , 
or about the mass of the Moon. The reason that this 
limit is not very restrictive is that the source stars are K 
and M-dwarfs, which are relatively small. Although the 
timescales of the events are unknown, the optical depth 
can still be estimated. In the limit of unresolved ( "spike" ) 
events, the optical depth becomes, 



Tobs 



spike 



1 



N*N d min(l,w|)' 



(Spike Events) (11) 



where Nsp^ is the number of spike events having A > 
1.34, Nd is the number of epochs per light curve, and Ut 
is, 
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(12) 



and At is the minimum detectable magnification. In the 
case of the six spike events, Ut > 1, and thus the optical 
depth is simply r obs = N spikc N^ 1 N ( [ 1 . Adopting N d = 25, 



r spikc = (2.9 ±1.2) x 10" 



(13) 



where the error is that due solely to Poisson statistics. 
The optical depth to the single resolved event can also 
be estimated using equation (0) , simply by setting Agp^e 
equal to the number of data points on the event light curve 
with magnification > 1.34, restricting attention to those 
taken in the F814W filter. Since there are six such points, 
this results in the same optical depth as T sp ike, in excel- 
lent agreement with the optical depth estimated from the 
usual formula. Thus the optical depth due to the spike 
events is comparable to that contributed by the one re- 
solved event. If the spatial and kinematic properties of 
the two samples are the same, this implies a similar mass 
in each component. 

4. BINARY LENSES, IONIZATION, AND EVAPORATION 

4.1. Binary Lenses 

All six of the spike events discovered by SM22 are well 
characterized by a flat light curve (to within the photomet- 
ric errors) with one deviant data point. As pointed out by 
SM22, the fact that the light curve exhibits no features 
other than the one deviant point implies that the planet 
must be quite distant from any parent star, or the lensing 
influence of the primary would be detectable. The mini- 
mum separation compatible with the data depends on the 
mass of the primary and the photometric accuracy, a p. 
Formally, the lensing behavior of the planet bound to a 
parent star is described by the formalism of binary lenses 
(Mao & Paczynski 1991; Gould & Loeb 1992). However, 
when the projected separation of the planet and star sep- 
aration in units of 6*e is much larger than the Einstein 
ring radius of the system, the magnification structure is 
well described by the superposition of two point masses 
(di Stefano & Mao 1996, di Stefano & Scalzo 1999b) sep- 
arated by a distance of d = (d 2 — l)/d, where d is the 
instantaneous angular separation of the planet and star in 
units of 6*e (Gaudi & Gould 1997). The detection proba- 
bility is roughly the probability that the source trajectory 
will pass within ut of the primary lens multiplied by the 
probability that the deviation due to the primary occurs 
during the observation window. Assuming that the mass 



of the secondary is much smaller than the mass of the 
primary, that the times of maximum magnification of the 
planetary events are uniformly distributed in the observa- 
tion window, and normalizing all distances to the mass of 
the primary, the detection probability is, 



Pdet(d) 




Au 



if d < ut, 
3in(Kt) if u T < d < u T + Ait, 

if d > ut + An, 

and Au = At/is, and is the timescale of the primary. In 
order to convert from the detection probability as a func- 
tion of d to the detection probability as a function of the 
semi-major axis a of the planet, it is necessary to convolve 
Pdct(d) with the probability of d given ap = a/ Re, 



Pdet{ap) 



dd P det {d)P{d; a P ) 



where 



P{d;a P ) = — I 1 
ap 
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(Gould & Loeb 1992). Pdet(ap) can be converted to phys- 
ical units by adopting a value for Re- Figure 1 shows 
-Pdct(a) as a function of a assuming a primary of mass 
M = 0.33M Q and A T = 1 + erp with a P = 1%, 5%, and 
10%. The 95% confidence level (c.l.) lower limit on a is 
the value where P dct {a) = 1 - 0.05 1 / 6 ~ 40%. I find 



a > 7.5 AU, 



(95% c.l.) 



(17) 



for ap = 10%. Figure 1 also shows the results of com- 
puting Pdct(ap) adopting the full binary formalism (Witt 
1990), which gives nearly identical results. 

The lower limit given in equation (|]j]) is somewhat 
model dependent; adopting a smaller mass for the primary 
would lower the limit. However, the estimate is conserva- 
tive, as I have assumed that the magnification from the 
primary must rise above At for it to be detected, when 
in fact, the cumulative effect of the curvature due the pri- 
mary on the light curve may make it detectable for sep- 
arations considerably larger than those calculated using 
equation (|l4|). 

Note that de la Fuente Marcos and de la Fuente Marcos 
(2001) give a lower limit on a that is a factor of ~ 2 smaller 
than equation ([17]) . This is due to the fact that they adopt 
a rather different - and considerably more conservative - 
detection criterion than assumed here. They require that 
d < ut for the primary to be detected. This requirement 
gives the absolute lower limit on d in the sense that the 
primary will be detectable in essentially all trajectories if 
d < ut- However, for larger d, a fraction P det (d) of trajec- 
tories should show evidence for a primary. Thus the lack 
of detections in six events can be used to place a more 
stringent lower limit. 

4.2. Ionization 

In dense stellar environments, such as the core of M22, 
planetary systems cannot be treated as isolated. The large 
number density and velocity dispersion imply that a star 
is likely, within the lifetime of the cluster, to encounter 
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another star. Such encounters will inevitably lead to strip- 
ping off of weakly bound outer planets, i.e. ionization. 

The average number of encounters a given planetary sys- 
tems has with a field star per unit time is given by l/t enc , 
where t cnc is the encounter time (Binney & Trcmainc 
1997), 



1 + 



GM* 
2a 2 a 



(18) 



where v is the local number density. For a = 1 AU, 
and assuming parameters appropriate for M22, t cnc = 

4.3 Gyr (v/lAl x 10 4 pc~ 3 )~\ Thus, in the core of M22, 
every star is likely to experience an encounter with another 
star with an impact parameter < 1 AU during the lifetime 
of the cluster. According to Heggie's Law, soft binaries 
will tend to get disrupted by such encounters. The binary 
is soft if a > a h / s , where 



2GM P 



J h/S 



3xlO _3 AU 



M P \ 
0.25Mj I 



11.4km s" 1 



( 19 ) 

and I have assumed the mass of the planet, M p , is much 
smaller than the field and parent star masses, which I have 
assumed to be equal. Thus essentially all planetary sys- 
tems are soft, and will be disrupted on a time scale t enc . I 
therefore write the ionization probability as a function of 
a as, 

P ion (a) = l-e- T °/*-H (20) 

This probability is shown in Figure 1 for three different 
values of the number density. Equations ( |l8| ) and (p(j|). 
although crude, agree well with more detailed and real- 
istic calculations (Smith & Bonnell 2001; Bonncll ct al. 
2001; Davies & Sigurdsson 2001). For the core of M22 
[y = 1.41 x 10 4 pc~ 3 ), the upper limit to the separation 
of a bound planet is, 



< 0.3 AU. 



(95% c.l.) 



(21) 



Comparing equations ( ^l|) and (|l7|), it is clear that the 
spike events cannot be due to planets bound to stars in 
the core of M22. 

There are several caveats. First is that I have not con- 
sidered the possibility of exchanges due to encounters, i.e. 
where the planet is transferred from the original parent 
star to the field star as a result of the encounter, rather 
than being ionized completely. The timescale for exchange 
is (Hut & Bahcall 1983; Hut 1983) 



3a 5 a 



20twG 3 M 3 ' 



(22) 



where I have again assumed that the mass of the planet is 
much smaller than that of the field and parent stars, with 
the latter two masses to be equal. For conditions in the 
core of M22, t cx = 1.1 Gyr(a/AU)(z//1.41 x 10 4 pc" 3 )- 1 . 
The ratio of the exchange to encounter timescales is given 

where I have neglected the term in t enc due to gravita- 
tional focusing. Note that t cx /t cnc is independent of v. 



(23) 



For a ^ 1.6 AU, the rate of encounters will exceed the 
rate of exchanges, and the net effect of encounters will be 
to ionize planets completely. Since this is safely below the 
lower limit on a set in §4.1, exchanges can be ignored. 

The second caveat is that I have not considered pro- 
cesses that may create planetary systems, such as three 
body interactions (Heggie 1975) and tidal capture (Fabian, 
Pringle, & Rees 1975, Press & Teukolsky 1977; Lee & Os- 
triker 1986). In fact, it can be shown that all such pro- 
cesses are extremely subdominant, and the net result of 
dynamical interactions on soft binaries is disruption (see 
Appendix 8.B of Binney & Tremaine 1997). 
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FIG. 1 The long-dashed curves show the ionization probabil- 
ity Plow as a function of semi-major axis assuming a cluster age 
of To = 12 Gyr, a velocity dispersion of a = 11.4 km s — 1 , and stellar 
densities of (left to right) v = 1.41 X 10 4 , 10 3 and 10 2 pc~ 3 . The 
solid, dotted, and short-dashed curves show the primary detection 
probability P^ ct as a function of the semi-major axis of the planet 
for photometric errors of op = 1%, 5%, and 10%, respectively. The 
triangles show the detection probability for <rp = 10% using the full 
binary-lens formalism. The bold solid curve is the probability of not 
being ionized and not being detected, i.e. (1 — Pi on )(l — ^det); f° r 
o-p = 10% and v = 10 2 pc~ 3 



4.3. Evaporation 

Some fraction of the planets that are ionized from their 
parent stars will have, upon ionization, speeds exceeding 
the local escape speed, which for the core is v e ~ 23 km s~ 4 
(Peterson & Cudworth 1994). These planets will rapidly 
escape the cluster on a timescale of order the crossing time, 
icross = T t /v e ~ 1 Myr. However, this will constitute a 
very small fraction of all ionized planets. The majority of 
the ionized planets will remain bound to the cluster with 
velocities similar to the velocities of their parent stars. 
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Over the lifetime of the cluster, these free-floating plan- 
ets will undergo many encounters with the cluster's stellar 
population. These encounters will tend to drive the system 
toward equipartition, so that cr 2 ~ (M*/M p )(7 2 , where er* 
is the velocity dispersion of the stars. If the cluster mass 
is dominated by the stellar component, i.e. the total mass 
in planets is negligible, then the potential will be deter- 
mined by the stars, and the local velocity dispersion a of 
the cluster will be cr* . In this case, it is clear that equipar- 
tition cannot actually be realized, since M*/M p > 4, and 
hence the planets will attain velocities larger than the lo- 
cal escape speed, a p > 2a = v e , and thus escape from the 
system. 

The timescale over which the planets get ejected from 
the system is the relaxation time (Spitzer 1987), 



irelax = 0.065 



vMlG 2 In A' 



(24) 



where In A is the Coulomb logarithm which I will take to 
be In A = ln(0.4iV,) = ln(47w r 3 /3) ~ 10. Here N* is the 
total number of stars in the system. For the core of M22, 
ireiax = 0.33 Gyr. Thus, assuming a constant density and 
velocity dispersion, the core has experienced ~ 35 relax- 
ation times during the lifetime of the cluster. Therefore all 
planets should have long since evaporated from the cluster 
core. 

Note that equation (M) is the relaxation time for the 
stars, which is appropriate for the planets only if the mass 
in planets is much smaller than that in stars. In the case 
when the two mass components are comparable, one might 
be tempted to adopt the equipartition timescale (Spitzer 
1969) 



«eq 



4V3nG 2 M*M p is In N 



'•relax- 



(25) 



However, as argued before, the system cannot actually 
achieve equipartition, because M*/M p > 4. Thus the es- 



timate m equation ( 
ration timescale wil 



breaks down, and the true evapo 
e considerably smaller. Regardless, 
if the mass in planets were comparable to that in stars to- 
day, M22 should exhibit an relatively large mass-to-light 
ratio, since planets of M p ~ 0.25Mj emit very little light. 
Peterson & Cudworth (1994) find a global mass-to-light 
ratio of M/L <; 1 for M22, which is anomalously low for 
globular clusters, and not compatible with a substantial 
population of planets by mass. Therefore equation (24) 
should be applicable. 

4.4. Planets in the Halo of M22? 

The results of §§4.1, 4.2 & 4.3 show that a significant 
population of planets, either bound or free-floating, can- 
not exist in the core of M22 (and explain the spike events). 
However, the relevant timescales, t Gnc and t rc i ax , are both 
inversely proportional to the stellar density v, which drops 
precipitously outside the cluster core, typically as v cx r -3 , 
where r is the radial distance from the center of the cluster. 
For the low density outskirts of the cluster, the encounter 
and relaxation timescales can be larger than the lifetime 
of the cluster, and bound and free floating planets can 
survive. 



Figure 1 shows t cnc for stellar densities of v = 10 3 pc~ 3 
and 10 2 pc~ 3 . For very low densities, bound planets are in 
principle compatible with the observed spike events. For 
example, for v — 10 2 pc -3 , there exists a range of a (albeit 
a small one) where bound planets could survive ionization 
and yet still satisfy the lower limit of a > 7.5 AU set by 
the lack of detection of the parent star. 

Similarly, for the outskirts of the cluster the relaxation 
time can exceed the lifetime of the cluster. Specifically, as- 
suming a is constant throughout the cluster, i re i ax > To for 
v <J 5.5 x 10 2 pc -3 . In these regions free-floating planets 
will evaporate very slowly. 
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FIG. 2 (a) The solid curve shows the physical density p(r) nor- 
malized to the central density po as a function of radius from the 
center of the cluster, r, normalized to the core radius r c , for a model 
with p(r) = p (l + r 2 /r;!)- 3 / 2 . For M22, p = 4.65 X 10 3 M Q pc~ 3 , 
r c = 1'.4, and the tidal radius is n = 29'.0 ~ 21r c (at the right edge 
of the figure). The dotted line shows the local velocity dispersion a 
normalized to the central velocity dispersions <ro. The short-dashed 
line shows the projected mass density E(r), where r now refers to 
the projected distance from the cluster center, (b) The long dashed 
line in the mass M(< r) interior to r, normalized to M(< rt), the 
total mass interior to the tidal radius, rt- The dot-short-dashed line 
shows /f ree , the fraction of free-floating planets that have survived 
evaporation. The dot-long-dashed line shows /bound: the fraction of 
systems that remain bound. 

Can the bound or free-floating planets in the halo of 
M22 explain the measured optical depth? The answer will 
depend on the competition between the declining contri- 
bution to the optical depth from the outer parts of the halo 
and the increasing longevity of bound and/or free-floating 
planets. Therefore it is necessary to assume a form for the 
density profile of M22. I will adopt the following form, 



p(r) =po[l + (r/r c ) 5 



-3/2 



(26) 
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Although a good match to observed cluster properties, this 
form has the disadvantage that the total enclosed mass 
diverges logarithmically. For practical purposes, I will de- 
note the total cluster mass as the mass interior to r t . From 
this form, the projected mass density, S(V), as function of 
the projected distance from the cluster center can be de- 
rived, 

Z(r) = 2p r c [l + (r/r c ) 2 ]~ 1 - (27) 

Figure 2 shows p(r), S(r), a(r), and M(r) as a function of 
r/r c . For M22, the concentration is c = log (r t /r c ) = 1.31, 
and r t — 20.4r c is just at the edge of the figure. 

From equation (ph, and given that size of the cluster is 
small compared to D s , the optical depth to microlensing 
is, 

r= 5—^(1 -a;), 28 

where x = Di/D s . Inserting values appropriate for M22, 
and adopting the form for given in equation (p7j), I 
find, 



Ttot 



10" 5 [1 + {r/r c y 



(29) 



Where r refers to the projected distance from the clus- 
ter center. Note that this overestimates the optical depth 
slightly, since the form of E(r) given in equation (^7|) is 
found from integrating the density along the line of sight 
from — oo to oo, when in fact the integral should be cut 
off at the tidal radius r t . The difference is negligible ex- 
cept near projected distances of r ~ rt, where the optical 
depth is extremely small anyway. The optical depth as a 
function r in shown in Figure 3. The HST observations 
of SM22 were concentrated in the inner ~ 2.5' of M22. In 
this region, t ~ 10~ 5 , which is a factor of ~ 3 times larger 
than the minimum optical depth estimated from the single 
resolved event (see §3.1). 

What fraction of optical depth can be contributed by 
planets? Assume some fraction f p of the cluster mass den- 
sity was originally in the form of planets. For simplicity, I 
will assume that this fraction is universal, i.e. the primor- 
dial mass density in planets is simply p p (r) = f p p(r). I 
will further assume that f p is sufficiently small that the dy- 
namics of the cluster is everywhere dominated by the stel- 
lar component, and furthermore that the cluster has not 
evolved significantly during its lifetime (almost certainly 
an oversimplification). The fraction of bound planets that 
will have survived ionization due to stellar encounters is 
simply 



/bound (r; a) = 1 - F lon (r; a) = e 



-Tb/te 



(30) 



where Pi on is given in equation (|2C|), and the dependence 
on r arises because t eric depends on v and a. An upper 
limit on /bound(?") is found by inserting a = 7.5AU, which 
is the lower limit on a set by the lack of detection of the 
primary in the spike events (§4.1). The fraction of free- 
floating planet that have survived evaporation from the 
cluster is, 

/frccM = l-^, (31) 
^rclax 

where £ e is the evaporation probability. Equation (|3l]) im- 
plies a linear dependence of the mass loss on time, as de- 
rived analytically for tidally truncated clusters by Spitzer 



(1987), and found numerically from the detailed evolu- 
tionary models (Chernoff & Weinberg 1990; Vesperini & 
Heggie 1997). For a relaxed system, £ e can be estimated as 
the fraction of stars in a Maxwellian velocity distribution 
that have velocities > v e — 2a, which is £ e = 7.4 x 10~ 3 . 
However, this is certainly an underestimate for the evapo- 
ration probability of the planets, because equipartition will 
drive them to significantly higher velocities than the mean 
a. I adopt the value of £ e = 0.156 for test masses given 
in Spitzer (1987), as derived from the models of Henon 
(1961). Note that this value of £ e is the global probability 
for the cluster referenced to the half-mass relaxation time, 
whereas for this calculation an estimate of the local evap- 
oration probability, referenced to the local relaxation time 
ireiax is necessary. I will simply assume the global value 
of £ e = 0.156 is appropriate, but note that this may be 
slightly in error. 
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FIG. 3 The solid line shows the total optical depth Ttot as func- 
tion of angular distance from the center of M22 in arcminutes for 
the model shown in Figure 2. The core radius (r c = 1'.4) and tidal 
radius (r t = 28'.9) are marked with downward-pointing arrows. The 
shaded box shows the inferred optical depth and error from the six 
spike events, T sp j ko = (2.9 ± 1.2) X 10~ 6 . The angular extent of the 
box is the longest width of the WFPC2 camera, 2'.5, and roughly 
corresponds to the region surveyed by SM22. The dotted line is the 
maximum optical depth from free-floating planets in such a model. 
The dashed line is the maximum optical depth from bound plan- 
ets. The dashed-dot line is the halo contribution to the total optical 
depth. 

The optical depth contributed from bound planets can 
now be estimated, 



4ttGDI _ 
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(32) 
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where the integral is along the line-of-sight over the extent 
of the cluster. The expression for Tf ree is the same, but with 
/bound replaced by / lree - The value for f p , the fraction of 
the original cluster mass in planets, is unknown. There- 
fore, I will simply consider the maximum optical depth, 
which is found by setting f p = 1. Note that this is truly 
a limit in the sense that it does not actually correspond 
to a self-consistent physical picture. I discuss this in more 
detail below. 

The resulting maximum optical depths to bound plan- 
ets, TboJ^, and to free-floating planets, r f ™^ x , are shown in 
Figure 3. For observations within the central ~ 2' of M22, 
r bou X „d = 6 - 4 x and T f ™f c x = 2.7 x 1CT 7 . These upper 

limits are ~ 5 and <~ 10 times smaller than the optical 
depth inferred from all six spike events, r sp ik c — 3 x 10~ 6 , 
and are only marginally consistent with the optical depth 
contributed by just one event. Note that such models re- 
quire that that outer 70% of the cluster be entirely com- 
posed of planets. 

How sensitive are the results to the adopted form for the 
density profile? For isolated clusters, the density should 
fall off as p(r) oc r~ 3 5 in the outer halo (Spitzer & Thuan 
1972), whereas tidally truncated clusters will exhibit some- 
what shallower profiles (Chernoff & Weinberg 1990). Both 
of these profiles are steeper than the p(r) oc r~ 3 behavior 
exhibited by the model adopted here. Indeed, the sur- 
face brightness profile of M22 falls off as p,(r) oc r~ 2 - 2 
in the outer halo (Trager et al. 1995), somewhat steeper 
than the X oc r~ 2 expectation of the adopted model as- 
suming mass traces light. In fact, the maximum optical 
depths for free-floating and bound planets are not very 
sensitive to the density profile. To first order, the frac- 
tions /bound and /f ree are inversely proportional to p(r). 
Therefore, to first order, the dependence of the maximum 
optical depth on the density structure cancels out, and 
only the dependence of the velocity dispersion a(r) on the 
radius enters into the calculation of the maximum opti- 
cal depth. Since this is a relatively insensitive function of 
the density structure, the maximum optical depths calcu- 
lated previously are reasonably robust to variations in the 
space of realistic density models. I have also computed 
the optical depth for a Plummer model which has a den- 
sity structure p{r) oc [1 + (r/r c ) 2 ]^ 5 / 2 . This yields a total 
central optical depth almost identical to that for the fidu- 
cial model, t ~ 10~ 5 . The maximum optical depths are 
r^ nd = 2.8 x 10~ 7 , and r™ a c x = 1.4 x 10~ 7 , about two 
times smaller than for the r~ 3 model, and 10 and 20 times 
smaller than inferred from the spike events. The true den- 
sity profile is likely bracketed by the Plummer and r~ 3 
models. 

It should be emphasized that the procedure of setting 
f p = 1 in the derivation of these maximum optical depths 
is not entirely self-consistent. The fractions /bound and 
particularly /f ree , were derived under the assumption that 
the dynamics are dominated by the stellar component. If 
this is not the case, then both the encounter time i C nc and 
relaxation time i rc iax may be considerably larger. Precise 
predictions of the behavior of the system in the presence of 
a substantial planetary component by mass would require 
a more sophisticated treatment than that presented here. 
However, based on the determination of the mass-to-light 
ratio of M22, it is certainly true that the core of M22 can- 
not be currently be dominated by planetary-mass bodies. 



Therefore, the conclusions of §§4.1-4.3 are secure: essen- 
tially regardless of its past dynamical evolution, a substan- 
tial population of planets cannot exist in the core of M22. 
As shown in Figure 3, within a projected radius of 1' from 
the cluster center, the halo contributes an optical depth 
of Thaio = 3 x 10~ 6 for the fiducial density model. This is 
also true for the Plummer model. Thus, reasonable mod- 
els can reproduce the observed optical depth only if the 
entire halo is composed of planets. This is essentially ex- 
cluded by measurements of the surface brightness (Trager 
et al. 1995) and mass-to-light ratio (Peterson & Cudworth 
1994) of M22. 

5. GALACTIC FREE-FLOATING PLANETS? 

The results of §4 strongly suggest that the spike events 
cannot be due to microlensing by lenses associated with 
M22. Therefore, if one were to continue with the inter- 
pretation that these events are due to microlensing, and 
assuming that the direction of M22 constitutes a generic 
line of sight toward the Galactic bulge, the lenses con- 
stitute a free-floating Galactic population. What do the 
observations of SM22 imply about such a population under 
this assumption? 

The optical depth inferred is independent of the location 
and nature of the lenses, therefore r sp ike = (2.9 ±1.2) x 
10~ 6 , as before. The upper limit on the mass of the lenses 
is now less certain, but assuming the lenses have typical 
bulge distances and kinematics, I find M p <; 0.5Mj - still 
firmly in the planetary regime. The lower limit on the 
separation from the lack of detection of the primary is 
a ;> 6.3 AU. Thus the planets must still have relatively 
wide orbits to escape detection if they are bound to par- 
ent stars. The MACHO collaboration has measured the 
optical depth toward the Galactic bulge based on two dif- 
ferent analysis methods. They find t = (2.0 ±0.4) x 10~ 6 
centered at (l,b) = (3°.9, — 3°.8) using clump giants as 
sources (Popowski et al. 2001), and r = (3.2 ±0.5) x 10~ 6 
centered at (2°.68, — 3°.35) using difference image analysis 
(Alcock et al. 2000). Note that these optical depths are 
mutually exclusive from that implied by the spike events, 
as the standard analysis techniques are not sensitive to 
events with timescales <J 1 day. The optical depth is 
a relatively strong function of position. The models of 
Han & Gould (1995) indicate that the optical depth at 
the position of M22 should be 3 — 4 times smaller than 
the values measured at the positions reported by MA- 
CHO. These estimates are necessarily model dependent, 
however, for definiteness, I will assume that the optical 
depth toward the MACHO fields and Baade's window, 
(I = l°,b = —4°), is three times smaller than toward 
M22. Thus the six spike events imply an additional op- 
tical depth of tbw = 3T sp ik c = 9 x 10~ 6 . Assuming that 
these planetary lenses have similar spatial and kinematic 
distributions as the bulge microlenses, this corresponds to 
> 1800(M p /0.5Mj)~ 1 planets per bulge microlens. 

These planets cannot be part of a halo population. The 
expected contribution to the optical depth toward M22 
from a standard, spherical, singular, isothermal halo is 
(e.g., Sackett & Gould 1993) 

^ = 5Xl °"(2^n7^) 2 (33) 
where is the asymptotic halo circular speed. Intro- 
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ducing a core to the halo would decrease this estimate. 
Even if only one of the spike events was due to microlens- 
ing by a halo object, this would imply an optical depth in 
planetary-mass objects toward the Large Magellenic Cloud 
in conflict with combined EROS and MACHO limits (Al- 
cock ct al. 1998). 

These planets also cannot be part of the thin or thick 
disk. The optical depth to an exponential disk is (Gould, 
Miralda-Escude & Bahcall 1994) 



be estimated as 



Tdisk = 



2ir GZ d D 2 s 
he 2 



where E<j is the local surface density of the disk, 

a = D s (\b\/h-l/R d ), 
g(a) = a- 2 [l - 2a" 1 + (1 + 2c*- 1 ) exp(-a) 



(34) 



(35) 
(36) 



and h and Rd are the scale height and scale length of the 
disk, respectively. I adopt Rd — 2.5 kpc (see Sackett 1997, 
and references therein). The optical depth is maximized 
for this value of Rd and the line of sight to M22 (b = — 7°.6) 
for h = 0.6 pc. The total observed surface density of mat- 
ter is 40 M Q pc- 2 (Gould, Bahcall, & Flynn 1996, Zheng 
et al. 2001), while the surface density of all matter be- 
tween ±1.1 kpc and the plane is 71 M Q pc -2 (Kuijken & 
Gilmorc 1991), leaving room for an additional 31 M Q pc -2 
of dark matter between ±1.1 kpc and the plane. This cor- 
responds to S,j = 37 M Q pc~ 2 , and thus the maximum 
optical depth toward M22 from a disklike component is 
Tdisk = 4.3 x 10~ 7 , about a factor of 7 smaller than that 
implied by the spike events. 

Finally, Binney, Bissantz, & Gerhard (2000) argue that 
the optical depth to resolved events measured by MACHO 
toward the Galactic bulge using difference image photom- 
etry, r = (3.2 ±0.5) x 10~ 6 , is already difficult to reconcile 
with our knowledge of Galactic structure (but see Sev- 
enster & Kalnajs 2001). Thus there is no room for the 
additional contribution of tbw ~ 9 x 10~ 6 , which would 
be required if the free-floating planets were smoothly dis- 
tributed. 

These arguments are, of course, indirect. However, there 
are several ways of unambiguously detecting or ruling out 
such a free-floating planet population. The (unpublished) 
results of the MACHO spike analysis toward the Galac- 
tic bulge would likely answer this question definitively. 
Performing spike analyses on the OGLE-I or OGLE-II 
databases (Udalski et al. 1993; Udalski, Kubiak, & Szy- 
manski 1997; Udalski et al. 2000) might prove more diffi- 
cult, due to the lack of contemporaneous color information. 
However, during the OGLE-III phase the sampling rates 
will be increased for some fields, enabling the resolution 
of considerably shorter time scale events. Finally, some 
events should be present in the databases of the follow-up 
microlensing collaborations. For example, the PLANET 
collaboration (Albrow et al. 1998) monitored N e <~ 60 
events toward the Galactic bulge with median sampling 
intervals of ~ 1 hour (Albrow et al. 2001, Gaudi et al. 
2002), sufficient to resolve essentially all events caused by 
planets capable of producing the M22 spike events. The 
number of expected events in the PLANET database due 
to the implied free-floating planet population can roughly 



A oxp = 



2A e (N*) r BW (At) 
7T (*b) 



(37) 



where (A*) ~ 10 3 is the average number additional stars 
on each frame, (At) ~ 40 days is the average duration 
that each event was monitored, and (ts) ^ 1 day is the 
average planetary event timescale. Thus A cxp ;> 15 events 
should be present in the PLANET database if the implied 
population of Galactic free-floating planets is real. 
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FIG. 4 (a) Differential distributions of the magnitude difference 
Am | for microlensing (dotted line) and a distribution that is uni- 
form over the observed range of |Am| (dashed curve). The curves 
are arbitrarily normalized. The vertical segments are the observed 
Am| for the six spike events, (b) The cumulative distributions of 
Am|, normalized to |Am| = 0.32 (or A = 1.34). The shaded his- 
togram is for the six observed events, the dotted histogram is the 
expected microlensing cumulative distribution, and the dashed his- 
togram is for the uniform distribution. The KS-statistic probabilities 
are indicated. 

All of these constraints can be avoided if the planets are 
not smoothly distributed. As suggested by de la Fuente 
Marcos & de la Fuente Marcos (2001), a dark cluster of 
planets along the line of sight toward M22 (but not associ- 
ated with M22) could reproduce the observed spike events 
without violating any of the limits on Galactic structure. 
The minimum mass of the dark cluster required to repro- 
duce the observed optical depth is, 



M dc 



7 x 10 4 M 







1 



(38) 



assuming the cluster radius subtends an angle > 2'. 5. 
Thus, unless the cluster is very close to us (Di <; lkpc), the 
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dark cluster must be quite massive, M^ c J> 10 M . See de 
la Fuente Marcos & de la Fuente Marcos (2001) for a more 
thorough discussion of the dark cluster interpretation. 

6. DISCUSSION 

The arguments of the previous sections lead to the con- 
clusion that the only explanation for the origin of the spike 
events that is consistent with all available observations 
and theory is a dark, massive structure composed of light 
M lMj compact objects coincidentally along our line 
of sight to M22. Such an explanation seems ad hoc at 
best. The simplest alternative is that these events are 
not due to microlcnsing. While it is not my intention to 
attempt to provide a definitive explanation for the spike 
events, in light of the evidence against the microlensing 
interpretation, it seems worthwhile to review and reassess 
the reasons why SM22 favored the microlensing interpre- 
tation over other sources of variability. 

The arguments in favor of the microlensing interpre- 
tation are that (1) the light curves are constant except 
for the spike event, (2) the colors and magnitudes of the 
sources are consistent with an unbiased sample of stars 
in the CMD and (3) the distribution of magnifications is 
consistent with that expected from microlensing. The first 
reason is compelling, and rejects a large fraction of known 
variables. However it is important to note that the SM22 
observations are restricted to 25 epochs over the course of 
100 days, and thus only exclude variables with duty cycles 
;> 4%. The fact that the distribution of source star mag- 
nitudes and colors is consistent with an unbiased sample 
of stars may not be very constraining, because the num- 
ber of proposed microlensing sources is quite small. Thus 
subtle biases are difficult detect. Also, it is possible, per- 
haps likely, that low-mass main sequence stars are those 
that are most likely to exhibit the required behavior: light 
curves that are constant to a few percent > 96% of the 
time with brief brightcnings of ~ 50% for < 4% of the 
time. If this were the case one would naturally expect the 
distribution of such sources to roughly trace the CMD. 
Finally the third argument, that the distribution of mag- 
nifications follows that expected from microlensing, is also 
not very constraining because of small number statistics. 
Figure 4 shows the cumulative distribution of | Am| for the 
six spike events. For microlensing, the cumulative distri- 
bution of magnifications normalized to A = 1.34 is simply 
P(> At) = where ut is the radius of the A = At 
magnification contour, and is given in equation ( [12] ) for 
point sources. The expected differential and cumulative 
distributions of |Am| for microlcnsing arc shown in Figure 
4. Also shown is a distribution which is uniform in \Am\ 
over the observed range. Both distributions describe the 
data equally well. 

At the distance to the Galactic bulge, and E(B — V) = 
0.33, these sources are brighter than the theoretical main 
sequence (Girardi et al. 2000) by about a magnitude. They 
would lie on the main sequence if at ~ 5 kpc. Alternatively 
the bulge stars could be significantly more reddened than 
assumed, E(B - V) = 0.73 versus E(B - V) = 0.33. Note 
that this would require a substantial amount of dust along 
our line of sight between M22 and the bulge. Without the 
original CMD, it is not possible to distinguish between 
these scenarios. However given that SM22 conclude that 



the number distribution of these source follows the dis- 
tribution of bulge sources, the latter scenario seems more 
likely. If these sources are in fact above the bulge main se- 
quence, it is possible that they sources are blends, e.g. M- 
dwarf/white dwarf binaries. Such a scenario may be able 
to account for the color and magnitude of these sources, 
and furthermore may explain the variability as well. 

7. CONCLUSIONS 

The primary conclusion of this study is that there is 
considerable circumstantial evidence that the six unre- 
solved ("spike") events detected by SM22 are not due to 
microlensing, and therefore SM22 have not detected a sub- 
stantial population of free-floating planets. The chain of 
logic is as follows: 

(1) Planetary-mass lenses bound to parent stars must be 

separated by J> 8 AU to explain the observations, 
else the influence of the parent lens would have been 
detected. 

(2) In the core of M22, the encounter timescale is t cnc w 

4 (a/AU)" 2 Gyr. Thus all planets with separations 
;> 0.6AU have been ionized by random stellar en- 
counters over the lifetime of the cluster, To = 12 Gyr. 

(3) The relaxation timescale in the core of M22 is t ro i ax = 

0.33 Gyr. Thus all free-floating planets have evapo- 
rated from the core. 

(4) For reasonable assumptions, the maximum possible 

optical depth to surviving planets in the halo of M22 
falls considerably short of the observed optical depth 
of r - 3 x 10~ 6 . 

(5) If smoothly distributed, the mass in free-floating 

Galactic planets required to produce the observed 
optical depth is extremely difficult to reconcile with 
current knowledge of Galactic structure. 

The only logical alternative is a dark cluster of planets 
with total mass M ;> 1O 4 M that happens to be along 
our line of sight to M22. This explanation seems ad hoc 
at best. 

Although there are no obvious alternative astrophysical 
candidates for these spike events, considering the weight 
of the arguments presented here, it would seem prudent 
to study existing data in order to fully characterize stellar 
variability at the relevant levels. Similar databases with 
higher temporal sampling, such as the HST time series 
photometry of the globular cluster 47 Tuc (Gilliland et al. 
2000), may be able to address this question definitively, 
without requiring additional resources. 
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